\o ^LV-VU ^taA Co» Cw4- fifPC , 

Ww:M o^i^. c*^. ^ 

tJ^4aar ^^dU. “ ^o ^ 
QSc^Aa &*- S*f. % ~Qd A t 


Propagation of cosmic rays: nuclear physics 
in cosmic-ray studies 


Igor V. Moskalenko* 1 , Andrew W. Strong** and Stepan G. Mashnik* 


* NASA/Goddard Space Flight Center, Code 661 , Greenbelt , MD 20771 
f Joint Center for Astrophysics/University of Maryland, Baltimore County , Baltimore ; MD 21250 
** Max-Planck-lnstitut fur extraterrestrische Physik, Postfach 1603 , D-85740 G arching, Germany 
x Los Alamos National Laboratory ; Lor Alamos , A/M 87544 


Abstract. The nuclei fraction in cosmic rays (CR) far exceeds the fraction of other CR species, such as antiprotons, electrons, 
and positrons. Thus the majority of infoimation obtained from CR studies is based on interpretation of isotopic abundances 
using CR propagation models where the nuclear data and isotopic production cross sections in p - and a-induced reactions are 
the key elements. This paper presents an introduction to the astrophysics of CR and diffuse y-rays and discusses some of the 
puzzles that have emerged recently due to more precise data and improved propagation models. Merging with cosmology and 
particle physics, astrophysics of CR has become a very dynamic field with a large potential of breakthrough and discoveries 
in the near future. Exploiting the data collected by the CR experiments to the fullest requires accurate nuclear cross sections. 


INTRODUCTION 

The origin of CR have been intriguing scientists since 
1912 when V. Hess carried out his famous balloon flight 
to measure the ionization rate in the upper atmosphere. 
The energy density of relativistic particles (CR) is ~1 
eV cm“ 3 and is comparable to that of the interstellar 
radiation and magnetic fields, and turbulent motions of 
the interstellar gas. This makes CR one of the essential 
factors determining the dynamics and processes in the 
interstellar medium (ISM). The observations of the Small 
Magellanic Cloud f 1] by the EGRET (Energetic Gamma 
Ray Experiment Telescope) on board of the Compton 
Gamma Ray Observatory (CGRO) have shown that the 
CR are a Galactic and not a “metagalactic” phenomenon. 
Observations of the Large Magellanic Cloud [2], in turn, 
have shown that y-ray flux is consistent with CR having 
a density comparable to that in our Galaxy. 

Major cosmic accelerators are supernova remnants 
(SNRs), with a fraction of CR coming from pulsars, com- 
pact objects in close binary systems, and stellar winds. 
Observations of X-ray and y-ray emission from these ob- 
jects reveal the presence of energetic electrons thus tes- 
tifying to efficient acceleration [3]. The total power of 
Galactic CR sources needed to sustain the observed CR 
density is estimated at 5 x 10 40 erg s"" 1 which implies the 
release of energy in the form of CR of ~5 x 10 49 erg per 
supernova (SN) if the SN rate is ~ 3 per century. 

Propagation in the ISM changes the initial composi- 
tion and spectra of CR species (Fig. 1). The destruction 


of primary nuclei via spallation gives rise to secondary 
nuclei and isotopes which are rare in nature, antiprotons, 
and pions (n ± , 7T°) that decay producing secondary e ±i s 
and y-rays. CR are “stored” in the Galaxy for tens of mil- 
lions years before escaping into the intergalactic space. 

Although much progress has been made since the di- 
rect measurements in space have become possible, the 
detailed information refers only to the environment near 
to the sun. The CR source composition and CR propa- 
gation history are imprinted in their isotopic abundances 
while diffuse y-rays and synchrotron emission from dif- 
ferent directions cairy clues to the proton and electron 
spectra in distant locations. These are the only pieces of 
the universal puzzle that we have and exploiting them re- 
quires extensive modeling. 



FIGURE 1. Basic processes in the ISM. 
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FIGURE 2. Spectrum ( E 2 x Flux) of diffuse y-rays from the 
inner Galaxy as measured by the EGRET. Curves indicate in- 
dividual components: 7i°-decay (NN), electron bremsstrahlung 
(EB), inverse Compton (IC), and isotropic diffuse emission 
(ID). Adapted from [4]. 

INDICATIONS OF NEW PHENOMENA? 

The puzzling excess in the EGRET data above 1 GeV 
(Fig. 2) relative to that expected [4, 5] has shown up in 
all models that are tuned to be consistent with local nu- 
cleon and electron spectra [6, 7]. The excess has shown 
up in all directions, not only in the Galactic plane. An 
apparent discrepancy between the radial gradient in the 
diffuse Galactic y- ray emissivity and the distribution of 
CR sources (SNRs) has worsened the problem [6]. 

Positron fraction e + /(all leptons) in CR as measured 
by HEAT [8] also exhibits an excess above —8 GeV 
(Fig. 3) compared to predictions of the diffusion model 
for secondary production [9]. 

Secondary antiprotons are produced in the same in- 
teractions of CR particles with interstellar gas as £ + ’s 
and diffuse y-rays. Recent p data with larger statistics 
[12] triggered a series of calculations of the secondary p 
flux in CR. The diffusive reacceleration models have cer- 
tain advantages compared to other propagation models: 
they naturally reproduce secondary/primary nuclei ratios 
in CR, have only three free parameters, and agree better 
with K-capture parent/daughter nuclei ratio. The detailed 
analysis shows, however, that the reacceleration models 
underproduce p’s by a factor of ~2 at 2 GeV [13] (Fig. 4) 
because matching the B/C ratio at all energies requires 



FIGURE 3. Positron s/(all leptons) ratio in CR compared to 
calculations in a leaky -box model [10] (solid) and GALPROP 
diffusion model [9] (dashes). Adapted from [11]. 
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FIGURE 4. Spectrum of secondary antiprotons in CR as 
calculated in reacceleration (solid lines) and optimized (dots) 
models. The upper curves - interstellar, lower curves - modu- 
lated (<I> = 550 MV) to compare with data. Adapted from [7]. 


the diffusion coefficient to be too large. 

If these excesses are not a simple artefact, they may be 
telling us about processes in the ISM, in the Local Bub- 
ble, or signaling exotic physics (e.g., WIMP annihilation, 
primordial black hole evaporation), but also may indicate 
a flaw in the current models. 


COSMIC RAYS AND DIFFUSE y-RAYS 

The modeling of CR diffusion in the Galaxy includes the 
solution of the transport equation with a given source 
distribution and boundary conditions (free escape into' 
the intergalactic space) for all CR species. The transport 
equation describes diffusion and energy losses and may 
also include [ 14] the convection by a hypothetical Galac- 
tic wind, distributed acceleration in the ISM due to the 
Fermi second-order mechanism (reaccleration), and non- 
linear wave-particle interactions. 

The study of transport of the CR nuclear component 
requires the consideration of nuclear spallation, radioac- 
tive decay, and ionization energy losses. Calculation of 
isotopic abundances involves hundreds of secondary nu- 
clei produced in CR interactions with interstellar gas. A 
thorough data base of isotopic production and fragmen- 
tation cross sections is thus a critical element of propaga- 
tion models that are constrained by the abundance mea- 
surements of isotopes, p’s, and e + ’s in CR. 

As solar activity changes with a period of 1 1 years so 
does the intensity of CR, but in opposite direction: with 
an inverse correlation. The “solar modulation” is a com- 
bination of effects of convection by the solar wind, diffu- 
sion, adiabatic cooling, drifts, diffusive acceleration and 
affects CR below 20 GeV/nucleon. The theory of solar 
modulation is far from being complete [15]. The Ulysses 
spacecraft first provided measurements of the solar wind 
and magnetic field outside the ecliptic helping us to un- 





FIGURE 5. The EGRET sky: diffuse y-ray emission [5]. 


derstanding the global aspects of modulation, while Pio- 
neer and the two Voyagers have explored the outer solar 
system. Recently there appear some indications that Voy- 
ager 1, currently at 13.3 x 10 9 km (88 AU) from the sun, 
may be close to the termination shock, the heliospheric 
boundary. If true, in a few years it will become the first 
spacecraft ever to reach interstellar space. 

The diffuse continuum emission is the dominant fea- 
ture of the y- ray sky (Fig. 5). It is evidence of CR pro- 
ton and electron interactions with gas and the interstel- 
lar radiation field (ISRF), and is created via 7T°-decay, 
inverse Compton, and bremsstrahlung. This emission in 
the range 50 keV - 50 GeV has been systematically stud- 
ied in by OSSE (Oriented Scintillation Spectrometer Ex- 
periment), COMFIEL (Imaging Compton Telescope), 
and EGRET on the CGRO and in earlier experiments 
[4, 5]. The GLAST (Gamma-ray Large Area Space Tele- 
scope) will improve the sensitivity for the diffuse emis- 
sion by a factor of 30. 

Increasingly accurate balloon-borne and spacecraft 
experiments are demanding propagation models with im- 
proved predictive capability. Incorporation of the realis- 
tic astrophysical input increases the chances of the model 
to approach reality and dictates that such a model should 
be numerical. Besides the nuclear data, such a model 
has to include the detailed 3-dimensional maps of the 
Galactic gas derived from radio and IR surveys, the Local 
Bubble and local SNRs, the spectrum of the ISRF, mag- 
netic fields, details of composition of interstellar dust, 
grains, as well as theoretical works on CR acceleration 
and transport in Galactic environments. The most ad- 
vanced model to date, GALFROP, is a three dimensional 
model [6, 13, 16]. It is widely used as a basis for many 
studies, such as search for dark matter signatures, ori- 
gin and evolution of elements, the spectrum and origin 
of Galactic and extragalactic diffuse y- ray emission, he- 
liospheric modulation. The model calculates CR propa- 
gation for nuclei (JH to ^Ni), p’s, e^’s, and computes 
y-rays and synchrotron emission in the same framework; 
it includes all relevant processes and reactions. It is an 
excellent tool to cross-test various hypotheses [7, 17]. 



FIGURE 6. Boron/carbon ratio in CR as calculated in the 
reacceleration model. Lower curve - interstellar ratio, upper 
curve - modulated (<t> = 450 MV) to compare with data. 
Adapted from [7]. 

NUCLEAR PHYSICS IN CR STUDIES 

The abundances of stable ( 3 Li, 4 Be, 5 B, 21 Sc, 22 Ti, 23 V) 
and radioactive secondaries (^Be, j 3 A1, ^Cl, ^Mn) in 
CR are used to derive the diffusion coefficient and the 
halo size [16, 18]. The derived source abundances of 
CR may provide some clues to mechanisms and sites 
of CR acceleration. However, the interpretation of CR 
data, e.g., the sharp peak in the secondary/primary nu- 
clei ratio (Fig. 6), is model dependent. The leaky-box 
model fits the secondary/primary ratio by allowing the 
path-length distribution vs. rigidity to vary. The diffusion 
models are more physical and explain the shape of the 
secondary/primary ratio in terms of diffusive reaccelera- 
tion (distributed energy gain) in the ISM, convection by 
the Galactic wind, or by the damping of the interstellar 
turbulence by CR on a small scale. 

The secondary/primaiy nuclei ratio is sensitive to the 
value of the diffusion coefficient and its energy depen- 
dence. A larger diffusion coefficient leads to a lower ratio 
since the primary nuclei escape faster from the Galaxy 
producing less secondaries and vice versa. The abun- 
dance of radioactive secondaries (e.g., ^Be^Be) is sen- 
sitive to the Galactic halo size, the Galactic volume filled 
with CR. The larger the halo the longer it takes for ra- 
dioactives to reach us thus decreasing the ratio ^Be^Be. 

Current CR experiments, such as Advanced Composi- 
tion Explorer (ACE), Ulysses, and Voyager, deliver ex- 
cellent quality spectral and isotopic data. Meanwhile, 
isotopic production cross sections have for a long time 
been the Achilles’ heel of CR propagation models, and 
now become a factor restricting further progress. Inter- 
pretation of CR data requires massive calculations of 
isotopic production involving p’s and a’s, however, the 
widely used semi-empirical systematics are frequently 
wrong by a significant factor [22, 23]; this is reflected in 
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FIGURE 7. Cross section for the reaction ™ f Si(/vr) 26 Al. 
Calculations: semi- empirical systematics [19, 20], CEM2K 
[21], evaluated [22]. Adapted from [21]. 


the value of propagation parameters. Figs. 7 and 8 illus- 
trate the effect of isotopic cross sections on the derivation 
of the halo size from radioactive secondaries. Using the 
semi-empirical systematics leads to the huge error bars 
where the upper limits are consistent with infinite halo 
size [24]. Using the evaluated cross sections dramatically 
reduces the error bars and gives a consistent value: 4-6 
kpc [22]. Unfortunately, the evaluated cross sections in 
the required energy range are mostly unavailable. 

K-capture isotopes in CR (e.g., 23 V, ^Cr) can serve as 
important energy markers and can be used to study the 
energy-dependent effects such as diffusive reacceleration 
in the ISM and heliospheric modulation [25, 26, 27]. 
Such nuclei usually decay via electron-capture and have 
a short lifetime in the medium. In CR they are stable or 
live longer as they are created bare by fragmentation of 
heavier nuclei while their j3 + -decay mode is suppressed. 
At low energies, their lifetime depends on the balance 
between the processes of the electron attachment from 
the ISM and stripping. The probability of attachment is 
strongly energy-dependent, increasing toward low ener- 



FIGURE 8. Determination of the Galactic halo size based on 
radioactive secondaries: using semi-empirical systematics [24] 
- dashes, light shaded area shows the range consistent with all 
ratios; using evaluated cross sections - solid, heavy shaded area 
shows the range consistent with all ratios. Adapted from [22]. 



FIGURE 9. 51 V/ 51 Cr ratio in CR as calculated with/without 
K-capture. Adapted from [26]. 


gies, while the probability of stripping is flat. This makes 
the abundances of K-capture isotopes in CR energy- 
dependent (Fig. 9). Without K-capture, the ratio 51 V/ 51 Cr 
in CR would be flat because both nuclei are secondary. 
The electron K-capture 51 Cr(EC) 51 V increases the ratio 
at low energies. Reacceleration (energy gain) increases it 
even further. On the contrary, solar modulation flattens 
the ratio (Fig. 10), which is very steep in the ISM. 

Study of the light nuclei in CR (Li-O) allows us to 
determine propagation parameters averaged over a larger 
Galactic region, but the local ISM is not necessarily the 
same and the local propagation parameters may signif- 
icantly differ. The best way to study the local ISM is 
to look at isotopes with shorter lifetimes (e.g., 14 C) and 
heavy nuclei since large fragmentation cross sections 
lead to a small “collection area.” 

The CR source composition is derived from direct CR 
data by correcting for the effects of propagation, spalla- 
tion, and solar modulation. The elements with low first- 
ionization potential (FIP) appear to be more abundant in 
CR sources relative to the high-FIP elements, when com- 
pared with the solar system material (Fig. 11). This might 
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FIGURE 10. 5 1 V/ 5 1 Cr ratio in CR as calculated for different 

levels of solar modulation. Interstellar ratio corresponds to 
<X> = 0 MV. Adapted from [27]. 
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FIGURE 11 - (Galactic CR sources)/(solar system) abun- 
dances of different elements vs. FIP. Adapted from [28]. 


imply that the source material for CR includes the atmo- 
spheres of stars with temperatures ~ 10 4 K [29]. A strong 
correlation between FTP and “volatility” (most of low- 
FTP elements are refractory while high-FIP elements are 
volatile) suggests that CR may also originate in the in- 
terstellar dust, pre-accelerated by shock waves [28, 30]. 
n Na, 31 Ga, 37 Rb and some other elements Z > 28 break 
this correlation. CR data tend to prefer volatility over FIP, 
but uncertainties in the derived source abundances (cross 
sections!) prevent an unambiguous solution. 

Isotopic peculiarities of CR composition are also im- 
portant. For example, the jyNe/^Ne enrichment might 
tell us that CR are produced in cores of superbubbles [3 1] 
created by multiple correlated SNe. 


SCENT OF NEW PHYSICS 

Galactic and extragalactic space presents a test range 
where nature runs its numerous experiments continu- 
ously for billions of years. This is an arena where all fun- 
damental forces perform in an exotic show involving yet- 
to-be-discovered particles, new elements, giant nuclei 
bound by gravitation - neutron stars, and singularities 
- black holes, and engineering the largest-scale grid of 
structures in the universe. CR and diffuse y- rays, there- 
fore, could contain signatures of exotic physics, however, 
conventional CR present an enormous background for 
tiny exotic signals. 

The growing number of experiments forces us to the 
conclusion that the universe is dominated by the dark 
matter (DM) and dark energy. A preferred candidate for 
non-baryonic DM is a weakly interacting massive par- 
ticle (WIMP). The WIMP is the lightest neutralino x° 

[32] , which arises in supersymmetric models of particle 
physics, or a Kaluza-Klein hypercharge B 1 gauge boson 

[33] . Annihilation of neutralinos creates a soup of par- 
ticles, which eventually decay to ordinary baryons and 



FIGURE 12. Longitude profile of diffuse y- ray emission 
[7]. Lines (top to bottom): total, 7r°-decay, inverse Compton, 
bremsstrahlung, extragalactic background. 


leptons. The DM particles in the halo or at the Galactic 
center [34] may thus be detectable via their annihilation 
products (e + , p, d, y-rays) in CR [35]. The approach is to 
scan the SUSY parameter space to find a candidate able 
to fill the excesses in diffuse y-rays, p’s, and e + ’s over 
the predictions of a conventional model (as discussed 
above). Preliminary results of the “global fit” to the e +, s, 
p’s, and diffuse y- ray data simultaneously look promis- 
ing [36]. In particular, the DM distribution with two-ring 
structure allows the EGRET y-ray data and the rotation 
curve of the Milky Way to be fitted, while the ring radii, 5 
and 14 kpc, surprisingly well coincide with the observed 
rings of cold H 2 gas and stars, respectively. 

In terms of conventional physics, the spatial fluctua- 
tion of CR intensity may also provide a feasible explana- 
tion. The CR p data can be used to derive the Galactic av- 
erage proton spectrum (Fig. 4, optimized model), while 
the electron spectrum is adjusted using the diffuse y-rays 
[7]. The model shows a good agreement with EGRET 
spectra of diffuse y-ray emission (<100 GeV) from dif- 
ferent sky regions (Fig. 12). The increased Galactic con- 
tribution to the diffuse emission reduces the estimate 
of the extragalactic y-ray background [7, 37]. The new 
extragalactic background [37] shows a positive curva- 
ture (Fig. 13), which is expected if the sources are un- 
resolved blazars or annihilations of the neutralino DM 
[38, 39, 40]. The discrepancy between the radial gradient 
in the diffuse Galactic y-ray emissivity and the distribu- 
tion of SNRs can be solved [41] if the ratio H 2 /CO in the 
ISM increases from the inner to the outer Galaxy. The 
latter is expected from the Galactic metallicity gradient 

The difficulty associated with p may also indicate new 
effects. The propagation of low-energy particles may be 
aligned to the magnetic field lines instead of isotropic 
diffusion [13]. Our local environment (the Local Bubble) 
may produce a fresh “unprocessed” nuclei component in 
CR at low energy [42]; the evidence for SN activity in 
the solar vicinity in the last few Myr supports this idea. 
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FIGURE 13. Spectrum of extragalactic diffuse y-ray emis- 
sion [37]. Adopted from [39]. 


CONCLUSION 

Several topics are expected to become the subject of in- 
tensive studies in the coming years. PAMELA (Payload 
for Antimatter Matter Exploration and Light-nuclei As- 
trophysics) is designed to measure p’s, e ±9 s, and iso- 
topes H-C over 0.1-300 GeV. Future Antarctic flights 
of a new BESS-Polar instrument (Balloon-borne Exper- 
iment with a Superconducting Spectrometer) will con- 
siderably increase the accuracy of data on p’s and light 
elements. The AMS (Alpha Magnetic Spectrometer) on- 
board the International Space Station will measure CR 
particles and nuclei Z<26 from GeV to TeV energies. 
This is complemented by measurements of heavier nuclei 
Z > 29 by (Super-)TIGER (Trans-Iron Galactic Element 
Recorder). The future GLAST mission will be capable 
of measuring y-rays in the range 20 MeV - 300 GeV; be- 
sides other goals, it should deliver a final proof of proton 
acceleration in SNRs - long awaited by the CR commu- 
nity. A breakthrough on SUSY and high-energy interac- 
tions should come with operation of the new CERN large 
hadronic collider, LHC. Not surprisingly, the success of 
the state-of-the-art CR experiments depends heavily on 
the quality of nuclear data and especially p- and a- 
induced reactions at intermediate energies from tens of 
MeV - few GeV. Challenging these new frontiers is thus 
impossible without involvement of the nuclear physics 
community. 

This work was supported in part by a NASA Astro- 
physics Theory Program grant, US DOE, and the CRDF 
Project MP2-3025. 


REFERENCES 

1. Sreekumar, P., et al., Phys. Rev Lett., 70, 127 (1993) 

2. Sreekumar, P, et al., ApJ, 400, L67 (1992) 

3. Koyama, K., et al., Nature , 378, 255 (1995); Allen, G.E., et 
al., ApJ, 558, 739 (2001) 


4. Hunter, S.D., et al., ApJ , 481, 205 (1997) 

5. Moskalenko, I.V., et al., “Diffuse Gamma Rays: 

Galactic and Extragalactic Diffuse Emission,” in The 
Multiwavelength Approach to Unidentified Gamma-Ray 
Sources , eds. K.S. Cheng and G.E. Romero, Kluwer, 
Dordrecht, 2004, in press (astro-ph/0402243) 

6. Strong, A.W., et al., ApJ , 537, 763 (2000) 

7. Strong, A.W., et al., ApJ, 613, 962 (2004) 

8. Barwick, S.W., et al., ApJ , 482, LI 91 (1997); DuVemois, 
M.A., et al., ApJ , 559, 296 (2001) 

9. Moskalenko, I.V., and Strong, A.W., ApJ, 493, 694 (1998) 

10. Protheroe, R.J., ApJ, 254, 391 (1982) 

11. Coutu, S., et al., Astropart . Phys., 11, 429 (1999) 

12. Orito, S., et al., Phys. Rev Lett., 84, 1078 (2000); Maeno, 
T., et al. , Astropart. Phys., 16, 121 (2001); Beach, A.S., et 
al., Phys. Rev Lett., 87, 271101 (2001); Boezio, M., et al., 
ApJ, 561, 787 (2001) 

13. Moskalenko, I.V., et al., ApJ, 565, 280 (2002) 

14. Zirakashvili, V.N., et al., A&A, 311, 113 (1996); Seo, E.S., 
and Ptuskin, V.S., ApJ, 431, 705 (1994); Ptuskin, V.S., et 
al., 28th Int. Cosmic Ray Conf., 1929 (2003) 

15. Potgieter, M.S., Spa. Sci. Rev, 83, 147 (1998) 

16. Strong, A.W., and Moskalenko, I.V., ApJ, 509, 212 (1998) 

17. Moskalenko, I.V., et al., A&A, 338, L75 (1998) 

18. Ptuskin, V.S., and Soutoul, A., A&A , 337, 859 (1998); 
Webber, W.R., and Soutoul, A., ApJ, 506, 335 (1998) 

19. Webber, W.R., et al., ApJS, 144, 153 (2003) 

20. Silberberg R., et al., ApJ, 501, 911 (1998) 

21. Mashnik, S.G., et al., Adv Spa. Res., 34, 1288 (2004) 

22. Moskalenko, I.V., et al., 27th Int. Cosmic Ray Conf., 1836 
( 2001 ) 

23. Yanasak, N.E., et al., ApJ, 563, 768 (2001); Moskalenko, 

I. V., and Mashnik, S.G., 28th Int. Cosmic Ray Conf, 1969 
(2003) 

24. Strong, A.W., and Moskalenko, I.V., Adv Spa. Res., 27, 
717 (2001) 

25. Soutoul, A., et al., A&A , 336, L61 (1998) 

26. Jones, F.C., et al., Adv Spa. Res., 27, 737 (2001) 

27. Niebur, S.M., et al., J. Geoph. Res., 108, A 10, 8033 (2003) 

28. Meyer, J.-P, et al., ApJ, 487, 182 (1997) 

29. Cass6, M., and Goret, P., ApJ, 221, 703 (1978) 

30. Epstein, R.I., Mon. Not. Roy. Astron. Soc., 193, 723 (1980) 

31. Higdon, J.C., and Lingenfelter, R.E., ApJ, 590, 822 (2003) 

32. Jungman, G., et al., Phys. Rep., 261, 195 (1996); 
Bergstrom, L., Rep. Progr. Phys., 63, 793 (2000) 

33. Cheng, H.-C., et al., Phys. Rev Lett., 89, 211301 (2002) 

34. Gunn, J.E., et al., ApJ, 223, 1015 (1978); Stecker, 

F.W., ApJ, 223, 1032 (1978); Zeldovich, Y.B., et al., 

Sov. J. Nucl. Phys., 31, 664 (1980); Silk, J., and Srednicki, 
M., Phys. Rev Lett., 53, 624 (1984); Gondolo, P, and Silk, 

J. , Phys. Rev Lett., 83, 1719 (1999); Cesarini, A., et al., 
Astropart. Phys., 21, 267 (2004) 

35. Bergstrdm, L., et al., ApJ, 526, 215 (1999); Baltz, E.A., 
and Bergstrom, L., Phys. Rev D, 67, 043516 (2003); 
Gondolo, P, et al., J. Cosmol. Astropart. Phys., 1, 8 (2004) 

36. de Boer, W., et al., hep-ph/03 12037; astro -ph/04 08 272 

37. Strong, A.W., et al., ApJ, 613, 956 (2004) 

38. Ullio, P, et al., Phys. Rev D, 66, 123502 (2002) 

39. Elsasser, D., and Mannheim, K., Phys. Rev Lett., in press 
(2004) (astro-ph/0405347) 

40. Salamon, M.H., and Stecker, F.W., ApJ, 493, 547 (1998) 

41. Strong, A.W., et al., A&A , 422, L47 (2004) 

42. Moskalenko, I.V., et al., ApJ, 586, 1050 (2003) 



